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Abstract 

C/^ ' Elliptic soliton solutions, i.e., a hierarchy of functions based on an elliptic seed solution, 

, are constructed using an elliptic Cauchy kernel, for integrable lattice equations of Kadomtsev- 

Petviashvili (KP) type. This comprises the lattice KP, modified KP (mKP) and Schwarzian 
KP (SKP) equations as well as Hirota's bilinear KP equation, and their successive continuum 
limits. The reduction to the elliptic soliton solutions of KdV type lattice equations is also 
discussed. 

> 

\D . 1 The lattice KP and Hirota equations 

m 

\ The study of the discrete versions of soliton systems, i.e., systems given by integrable par- 

tial difference equations, have become in recent years a focus of attention in the theory of 
integrable systems. Among those systems, the discrete analogue of Kadomtsev-Petviashvili 
(KP) equations which define in three dimensional lattice, seem to form a universal class of 
systems. The first equation of this type was found by Hirota in [T^] and was referred to as 
DAGTE (Discrete analogue of generalised Toda equation) which is the bilinear equation 



(aie^i +a2e^^ +a3e^^)r-T = 0, (1.1) 

where the Hirota operators Di produce finite forward-and backward shifts, when acting on 
a pair of functions, in the corresponding lattice direction, i.e., 

e''V-5 = /(m + l,....)5(ni-l,--) ■ 

Special reductions of this equation, are obtained when the coefficients ai, 02, 03 satisfy the 
condition ai + 02 + as = 0, but the full equation is integrable in the sense of multidimensional 
consistency for arbitrary parameters, see [5]. Miwa [TH] reparametrized the equation in that 
restricted case, and hence in that form it is often referred to as Hirota-Miwa equatioi^. In 
this paper we will investigate a class of solutions of this and related three-dimensional lattice 
equations, comprising the following equations: 



^Sometimes the full equation (representing the blinear discrete KP equation) is also (in our view erroneously) 
referred to as the Hirota-Miwa equation. In fact, in [18] only the restricted case was considered, and generalized 
to a four-term equation which is nowadays referred to as the Miwa equation. 
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The bilinear lattice KP equation 

ff + ff-ff = Q- (1.2) 

The lattice KP equation 

{w — w){w — w) = (w — w){w — Hj) . (1-3) 
The lattice modified KP equation 

V — V V — V V — V 



V V 



(1.4) 



The asymmetric modified KP 



= 0. (1.5) 



The lattice Schwarzian KP 



{z — z){z — z){z — z) 
{z — z){z — z){z — z) 



(1.6) 



The notation in (ll.2p - (jl.6p is as follows: all dependent variables are functions defined on 
the multidimensional lattice with discrete coordinates {n,m,h) e Z'^ , e.g., / = f{n,m,h), 
and the elementary shifts in the three discrete directions are denoted by / = /(n + 1, to, /i), 
7 = f{n,m + l,h) and / = f{n,m,h + 1). The first equation ([Ll]) is Hirota's DAGTE 
after a change of independent variables and a point transformation. The other KP type 
lattice equations (|1.3p - (|1.6l) where established in [52] in a form containing parameters but 
which are equivalent to the normalised forms given above by point transformations. In 
particular, the Schwarzian KP (jl.6p was first given in this normalised form in [10] . and 
was subsequently recovered in [5] , whilst its geometric significance was explored in [Ul [TU 
[^. All equations (|1.3p - (ll.6p appear as distinct parameter choices of a generalized lattice 
KP equation given in ^22j, whilst the (|1.2[) can be viewed as a potential equation for the 
asymmetric lattice MKP equation. With regard to terminology, we prefer to reserve the name 
lattice (potential) KP equation for (|1.3|) rather than for (|1.2|) or for (|l.ip (as is commonly 
done in the literature) , because it is actually the former equation that is more directly related 
to the actual continuous KP equation whilst the latter two equations are more related to the 
bilinear form of the KP equation. Similarly, the equations (|1.4p and (|1.6p can be shown to be 
related directly to the continuous (potential) modified and the Schwarzian KP respectively 
equations through continuum limits, cf. [HI [31] ^T\, hence their names. It was recently 
established in '5^ that this list of five canonical equations exhausts all possible cases of 
octahedral type lattice equations, up to equivalence, which are multidimensionally consistent 
when embedded in the four dimensional lattice. 

Some of these equations have been quite widely studied. In particular, eq. (jl.l[) has 
arisen in many contexts over the last few years, notably in connection with affine Weyl 
group description of discrete Painleve equations, [25], and in connection with the represen- 
tation theory of Kac-Moody algebras [13]. An interesting connection was found with the 
spectrum of Bethe Ansatz states of quantum solvable models, the eigenvalues of which were 
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shown to obey specific versions of Hirota's bilinear equation, cf. [171 US- Furtliermore, 
related to this is a connection found in |24) between, on the one hand, pole-type (in the 
sense of Airault, McKean and Moser [T]) and soliton solutions of the lattice KP equation 
and a class of discrete-time many-body systems which effectively are gives rise to integrable 
correspondences (i.e., multi- valued dynamical maps) which are time-discretizations of the 
Ruijsenaars model [5^1 [22. On the other hand the equations describing these same classical 
discrete dynamical systems coincide exactly with the Bethe-Ansatz equations for the exci- 
tations of quantum integrable models solvable by the quantum inverse scattering transform. 
It seems from all these intriguing connections that the lattice KP equations exhibit a certain 
universality in their connection with various types of integrable models involving partial dif- 
ference equations, classical many-body systems and quantum solvable systems, a full picture 
of which is not yet understood at this juncture. 

In this paper, we consider a novel class of solution of these lattice equations namely 
soliton type solutions based on an elliptic seed solution. The construction exploits elliptic 
N X N' Cauchy-matrices, and involve also N' x N coefficient matrices, which in principle 
allows for a classification of the solutions according to the Schubert decomposition of a 
corresponding Grassmannian. Thus, our main result, which provides a general expression 
for what we accordingly call the elliptic {N, A^')-soliton solutions, could form a basis for a 
similar analysis as was performed in the case of the continuous KP equation by Kodama 
and Chakravarty, |16] , to describe the soliton taxonomy of the KP equation. Furthermore, 
our approach yields the various Miura type relations between the different equations and 
allows to study the various intermediate continuum limits, as well as reductions to lower 
dimensional equations, in a systematic way. 



2 Cauchy matrix scheme 

We will develop now a scheme along the lines of the paper [ini HO] where Atkinson and one of 
the present authors developed the elliptic soliton solutions for the ABS (Adler/Bobenko/Suris, 
cf. [3]) list of two-dimensional lattice equations, based on elliptic Cauchy matrices. That 
construction provided solutions for all equations of the ABS list up to Q3, whereas the case 
of Q4 was treated separately using a different approach, cf. [5] . 

In this Section we derive the basic relations, and in the next Section we will use these 
relations to find a general elliptic (iV, iV')-soliton solution for KP. 

2.1 Basic ingredients 

At this point we find it useful to introduce the Lame function 

*4(k) =: $4(K)e-^(«)" , (2.1) 

where 

cr(Ocr(K) 

with a is the sigma function and we have included an exponential factor together with the 
^-function, breaking the symmetry between the argument of the function and the suffix. 
Although most of the results of this paper can be obtained in terms of the ^-function alon^l, 
i.e. without adopting the exponential factor, the exponential factor does provide a certain 
regularisation of the functions involved which may affect their analytic behaviour. The basic 
identities for the ^ function are the following: 

^^{k)^s{X) = e'««*5+5(K)*5(A-«;)+e''^^^'4(At-A)*4+5(A) , (2.3a) 

= e'^^-^^+s{n)[ao+m + a^)-c{c+s+^)] , (2.3b) 
vi/^(k)*4(a) = ^d^+x) m)+a^) + cw-a^+^+^)] . (2.3c) 



^The inclusion of the exponential factor amounts to a specific gauge transformation on the quantities defined 
later on in the construction, and hence could be removed without affecting the main results. 
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in which we have introduced 

.. = ..(0 = c(e + ^)-c(e)-cw4^g|^. (2.4) 

Furthermore, we have the symmetry: k) — — . 

The starting point for our construction is the "bare" non- autonomous Cauchy matrix 

M° = (AfO^.).,,.i,...,jv ' + , (2.5) 

depending on a variable ^ depending hnearly on independent variables n, m, namely ^ = 
5o + JT-^ + ?7i£ + , with 5, e and A being the corresponding lattice 'parameters. We will 
assume that the set rapidity parameters {k;, i = 1, . . . , iV} is such that Hi + k^- / (modulo 
the period lattice of the cr- funct ion ) . 
Furthermore, defining 

= *5(k) , (Ik = *e(K), = *a(k), (2.6) 

and setting = ±K^-, we can derive from the basic addition formula (I2.3ap the 

following dynamical properties of the elliptic Cauchy matrix: 

MI^Vk'^ - *e(K, + At;-)*5(«j) 

= e''^(^'+^^)vl;^+i(«;, + 4)vI/,-(-AC,) + e''^'=^^P^+5(4)vI/5(K,) 
= M5p_«,e/«(-'+^^) + e''^-^vI/^+5(K;.)*?(«.) . 
We introduce now the plane- wave factors (i.e., discrete exponential functions) 

K'^) = (e-';(^)>_«) " (e-^(^)«<z_,) " (e-«^)«/_,) e^(«)>o,o.o('*), (2.7) 
K'*') = (e'^'^'^'^K') "™ (e"^)'^'/^') "^^^^''Vo,^^^ (2.8) 

and for the specific values k = k^, pi :— p{ni) and k' = nip Vj :— p{k'j) obeying the shift 
relations ~ ^ _ 

^^e^^-'p^^. , ^^e^^^'g-^, , ^ = e''-'='Z_,,, (2.9) 
Pi Pi Pi 

^=e''^«HPK')"' , !^^e^i^^'^iq,,)-' , ^^e-'^^H^^O^', (2-10) 

where the superscripts " ~ " ^ " and " ~ " denote the lattice shifts related to the 
shift by one unit in the variables n, m and h respectively, making use also of the relations 
^ = ^ + (5 , ^ = ^ + e and ^ = ^ + A which sit inside the coefficients rjs, rj^ and rjx- 

Now we can introduce the N- and iV'-component vectors 

r=iPi^d^^)),^,_^ , s={iy,^^iK'^))^^^^^^, (2.11) 

where pi = Pn,m,h{ni) and Vj ~ Vn.m,h{n'j) vly terms of which we can define now the 
"dressed" Cauchy matrix: 

M = (M,,,),,^i,...,^ , M,- , = PiMljV, . (2.12) 

As a consequence of the relations given earlier, and employing the definitions of the plane- 
wave factors, we can now describe the discrete dynamics as follows: 

Lemma 2.0.1. The dressed Cauchy matrix M, as defined in (|2.12p . obeys the following 
linear relations under elementary shifts of the independent variables n 

M ^ M ~rt^ , (2.13) 
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and under shifts of the variable m the similar relations: 

M = M -rs^, (2.14) 
and under shifts of the variable h the similar relations: 

M = M -rs'^, (2.15) 

where, as before, M , M and M denote the shifted Cauchy matrices. 

In what follows we will employ the relations (j2.13p , (|2.14p and (j2.15l) to obtain nonlinear 
shift relations for specific objects defined in terms of the Cauchy matrix M. 

2.2 The r-function and related basic objects 

Introduce now the T-function : 



T ^Tnmh= det (1 + MC) = det (1 + CM) , 

' ' NxN N'xN' 



(2.16) 



where, since M is in general not a square matrix, the N' xN constant matrix C is introduced 
to compensate for the discrepancy. The matrix 1 is either the N x N or N' x N' unit 
matrices respectively in both determinants. The latter identity is a consequence of the 
general Weinstein-Aronszajn formula: 



det 

NxN 



N' 

E 

1=1 



mici 



det (l 

N'xN' ^ 



Ik) , 



(2.17) 



where mj = {Mi^i)i=i n and c* — (Ci,/)j=i,...,Ar are the N' column-resp. row vectors from 

the matrices M and C. 

From the definition of the r-function and the relations for the dressed Cauchy matrix M 
we obtain: 



1 + Mc) ^ det il + (M -rs^)c\ 

J NxN I J 

{(1 + MC) 1- {1 + MCy^rs'^C ^ 



det 

NxN 

det 

NxN 

T det ( [l - (1 + MCy^rs'^c] ] 

NxN II J J 



Introducing 

= XcAO ^ C(a) + m + C(C) - C(C + a + /3) 

and using the fact that 

^~(K') 



(2.18) 



(2.19) 



where pfc is the diagonal matrix with entries p^' and K' = dia.g{K[, k'^). Here and 
in what follows the notation xs,k' denotes the diagonal matrix with entries xs,k'.j [j = 
1,...,N'). 

Now we have ^ 

(2.20) 



- = l-s^ [X5,K'r' C[l + MC]-^ r = Ws. 



The reverse fraction of Eq. (|2.20p can be computed by processing the same computation 

det^ (1 + MC) = det^ |l + (M + rs'^)c| 



NxN 

= det 

NxN 



NxN 



T det 

NxN 



{(1 + MC) l + (l + MC)-Vs^C } 
{ l + (l + MC)-Vs^C }. 
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Using also the fact that 



r = e" p-K .-r. / r^x r ^- [x-&,k\ r, 



where K = diag(Ki, kn) and then we have 



— = 1 — S 
T 



l + MC 



For arbitrary parameter a we introduce 

= l-s^[l + CM]-^C[xa,if]"'r, 



where we conchide that 



(2.21) 

(2.22) 

(2.23) 
(2.24) 

(2.25) 



We will now proceed the derivation of the dynamical relations of Ws and Vs. 
2.3 Basic linear relations 

In order to derive relations for the objects Va and Wa we need to introduce the iV- and 
A^'— component column- rcsp. row vectors: 



w„ = {l + MC)-\xc.,Ky^r , 
V = s'^{xp.K')-\l + CM)-^ . 

Using these equations, we can write the functions V and W in the forms 

Va = I — S^CUa , 
Wfi = 1- %Cr . 

Performing the following calculation: 



(2.26a) 
(2.26b) 



(2.27a) 
(2.27b) 



Ua = {l + MC)-'ixa,K)-'r = (1 + MCy'ixa.Kr'e'-' 



(1 + MC)uc 



{l + MC)-rsC 



CiK) + C(e + «) - C{S) -CiK + C + a) 



CiK) + + a+a) - CiK + ^ + a) 

c(+a) + c(e)-c(e + «) + c(^) \ 



\^ CiK) + + a+a) ^C{K + i + a)) 
Multiplying both sides by (1 + MC)^^ and introducing the vector 

Mo = (1 + MC)-V , 

we get the relation 



(2.28) 



Ua = -MqVq + Xa,5Ua ■ (2.29) 

A similar set of relations can be derived for the adjoint vectors (|2.27b|) which involves the 



adjoint vector to (I2.28p . namely 



*Mo = s'^{l + CM)-^ 



(2.30) 



and obviously these relations all have their counterparts involving the other lattice shift 
related to shifts in the discrete independent variables m and h instead of n. 

Summarising the results of these derivations, we have the following statement: 
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Lemma 2.0.2. The N- and N' -component vectors given in (j2.26l) . together with the ones 
given in (j2.28p and (j2.30p obey the following set of linear shift relations 

Ua = -UoVa + Xa,SUa , (2.31a) 

V = Wp *Mo - X0.-5 (2-31b) 
and a similar set of relations involving the shifts in the variable m is given 

Ua = -UoVa + Xa,eUa , (2.32a) 
= Wp *Wo - X0,-e (2.32b) 

and a similar set of relations involving the shifts in the variable h is given 

Ua = -UoVa + Xa,\Ua , (2.33a) 

= Wp *Wo - Xp,-x (2.33b) 
2.4 Basic nonlinear relations 

The nonlinear relations can be obtained by multiply s^C on the left hand side of (|2.31ap . 
We have 

S^CUa ~ Cu^Va + Xa,8S^ CUa 

s'^X5,K'CUa = -WoVa + Xa,s{l ~ Va) , (2.34) 

where wq = s^Cuq. The same relations can be obtained in the same way for other shift 
directions 

's^Xe,K'CUa = -W^Va + XaA'^ - ^a) , (2.35) 
S^X\,K'CUa = ~WqV a + Xa,\{'^ - Va) . (2.36) 

Combining (|2.34p and (|2.35p . we obtain 

V a{w - W) ^ Xa,5Va ~ Xa,EVa , (2.37) 

where w ~ C(0 ^ ~ n({S) — mC,{e) — h(^{j) . 

Similarly, multiplying the right hand side of (j2.31bP with Cr, we obtain 

Wp{w -w)= X0,-5Wp - X0.~eWp . (2.38) 
We now can write w — w in the following 

= ^ ( ^imtfiwp - Mm^wp] . (2.40) 

Wp \ ^p{i-s) ^p{(-5) 'J 

Introducing new variables Vq = ^aiO^a and = ^p{£^)Wp, we now have 

(VU-W) = ^ (paVa - qaVa) - [p-P^^P " 9-/3 W^) ■ (2.41) 

This equation has its counterpart involving the other lattice directions namely 

(W-W) = ^ (laVa ~ QaVa) = ^ [l-P^^p - Q-pV^p) • (2.42a) 
(w-w) = ^ [paVa - LVa) = ^ {p-P^^P " l-P^^p) ■ (2.42b) 
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Using (1231]), ((2:423 and (11329, we have 

Vc % 

or equivalently 



= , (2.43) 



From (|2.42a|) . if wc take a = — e we also have 



T T 



(2.44) 



(2.45) 



W/3 W/3 ~ ' 

These two equations are the "modified lattice KP" . 

Taking a = —S, (12.43^ becomes 
Taking /3 = (5, (piU)) becomes 

(|2.45p and (|2.46p are the ''asymmetric modified KP^\ 
Furthermore, for a = —S, we find that 

(w-w) = -q-sV-s/V^s, (2.47) 
{w-w) = -l-sV-s/^-s- (2.48) 

The combination of these two equations yields 

{w — w){w — w) = (w — w)(w — w) , (2.49) 

which is the ''lattice KP equation^' . 

Using (12251), we can write (12^47)) and in the forms 

(w-w) = -q-sl^, (2.50) 

T T 

(w-w) = ~l-sL^- (2.51) 
r r 



(lI;-W) = -q-e--. (2.52) 
T T 

The combination of (P^TI) and (1^3^ gives 

Z_5rf + g_efT - g_5rf = , (2.53) 

which is actually the "Hirota's DAGTE" . In the rational case, the summation of coefficients 
would add up to be zero. In the elliptic case, this condition is no longer to be the case. 
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2.5 Schwarzian KP variables 

From the previous relations for Ua we can now proceed as follows 

■ — ^ — ^ g _ _ 2J, 

" gK') + g(3) + go-gc + K' + i3) " 
~T (, ^ g5)-gp) + g^ + P)-gi + 6) \ ~ ,^ .7 . ^ ^ 

= s 1 + — — ; Mq = (1 - Vq) - Xl3.-S^i3.a 

\ c{K') + g(3) + go-g^ + K' + p)J > p. p. 

= Xa.sSp^a - (1 - Wp)Va , 

from which we get the following relation; 

WisVa = 1- X0-sSp,a - Xa,sSi3,a , (2-54) 

where 

S0,a = s^[xp,K'r'C[l + MC]-^[xa.Kr^r , (2.55) 

= *upC[xa^Kr^ r = s'^ [Xf3^K'r^ CUa . (2.56) 

Similarly we have 

W/jVa = 1 - Xl3~eSl3,a " Xa.sS/^.a , (2.57a) 

WpVa = 1 - X/3 -A'S'/J^Q - Xa.xSp^a , (2.57b) 

Note that if we multiply C [xo.k] ^ r from the right hand side of Eq. (|2.33b[) the computa- 
tion leads to Eq. (g^. 

Multiplying $0(0*5' ,9(0 to the left hand side of (|2.54p . we have 

+<^^+p[0^p{-5) + $,+^(0<i>a(<5) 

+ '^a+p{C){jPa ~ PaXa.liSp^a) ■ (2.58) 

Introducing the new variable 

^/3,a = $^,a(0(l-Xa./3'5/3,a) > (2.59a) 

we can write p. 581) in the form 

VqW^ = PaZp^a - P-pZfj^a ■ (2.60) 



Similarly we have 



Using the identity 



VqW;3 = QaZp^a ~ Q-fiZfj^a , (2-61) 

V^Wp = l^Zp^^ - l-fiZp^^ . (2.62) 



(VgWff) ^ (Vc^Wff) (V„W^) 

(Iw^) (?„W^) (V„W^) ' 

we can derive the equation 

[qaZp.a - q-pZp^a) [LZ p ,a - I -0 Z j3 ^a) {qaZ p^a - Z p ^a) 



(2.63) 



(2.64) 
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which is the ^^Schwarzian lattice KP equation" , first given the exphcit form in |10j . 

Remark: Let d denote the derivative w.r.t any independent variable on which solely the pi 
and all the Vj depend, but not any of the other ingredients in the elliptic soliton solutions. 
Hence d only acts on the pi and Vj . Assuming that the latter can be written in the following 
form: 

r = Re , = e^S where = (1, 1, 1). 
we can use the following relation for the Cauchy matrix 

dM = {dR)R'^M + MS^\dS). 

to derive the following expression for the action of d on the main variable Sa,i3 

dSc,,0 = [S^^aSC + CdRR^] u^. (2.65) 

This expression, in fact, is the analogue of the square eigenfunction expansion (in the sense 
of the seminal paper by Deift, Lund and Trubowitz of 1980, cf. ^ ) of the elliptic soliton 
solution of a continuous KP equation, e.g. by choosing d to represent the partial derivative 
w.r.t. one of the continuous independent variables of the KP equation. 

In the discrete case one can derive a similar equation to (j2.65p . Let" ~ " denote here the 
shift w.r.t any discrete variable on which solely the pi and Uj depend. Using the ingredients 
of the previous case, one can derive the formula: 



SaR — SaM — 



CRR^'-SS^'C M„. (2.66) 



which constitutes the discrete analogue to (j2.65p . Note that the shift " ~ " may also represent 
a composite shift, or a combined shift w.r.t. multiple variables, the derivation of (|2.66p only 
uses the fact that the shift distributes over products of functions of the discrete variable by 
f9 = f9 that the discrete variable enters the solutions via the pi and the Vj. There is, 
however, an important difference between the formulae (j2.65p and (j2.66p . namely that the 
left hand side of the latter involves also shifts of the eigenfunctions. We expect that these 
formulae may prove useful in the derivation of conservation laws for the discrete equations 
considered in this paper. 



2.6 Hirota form of the elliptic A^-soliton solution 

In this Section, we will derive some explicit formulae for the soliton solutions in terms of the 
T-function, which allow us to study their properties. Using the fact that the matrix M is 
actually a Cauchy matrix, the r-function can be explicitly computed by using the expansion 

N 

T = det (1 + MC) = 1 + ^ \B^.^\ + J2 

i—1 i<j 

where B = MC. 

Lemma 2.0.3. Cauchy-Binet formula: For an arbitrary N x M matrix A and M x N 
matrix B we have the following formula for the N x N determinant of the product ill}/ : 

( if M <N 

det {AB) = I det{A)det{B) if M^N 

^""^ I Ei<;i<...<iiv<A/det (A(i,..,A,)(i^,...,,„)) det (B(,^,...,,„)(i,..,Ar)) if M > N 

in which A.(x,..,Ar)(/j^_. denotes the matrix obtained by selecting the h, ...jIn columns from 
the matrix A and (^^^(i ^y) the matrix obtained by selecting the h, In rows from 

the matrix B. o 



BiA B 



B 



3,1 



B 



1,3 
3,3 



det(B) 
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Using Cauchy-Binet formula, we may express the r— function in the form 

N / N' 



)N I N' 
+ E E 
i<3 \h<k 

Ma, Mu, Mu, 
MjW Mji, Mji, 
Msi, M,i, Msi, 



Ma, Ma, 



Ci^i Ci^j 



N 



N' 



E E 

i<j<S \h<l2<h 



Ci,i Ciij Ci,s 

Chi Cl^j C12S 

Ci^i Ci^j Ci^s 



^ det M(i^2,...,N)(h,/2....,!iv) det C(;j,,,^...^;„)(i_2,, 



,7V)- 



(2.67) 



ll<l2<---<h 



where are the entries of the matrix C and M (i,2,...,n)(Ii,12,---,In) denotes the matrix ob- 
tained by selecting the (Zi, /2, ^at) columns from the matrix M and (^((^ j^,. ..,/«)(!, 2,..., w) 
is the matrix obtained by selecting the (^i, I2, In) rows from C. 

Using the Frobcnius formula for the relevant elliptic Cauchy determinants. Thus, from 



det 



+ K.j)a{n[ + k'j) 



(2.68) 



Introducing the notations 



e""'"' = Cr(Ki — Kj) 



the Hirota formula for the r-function thus takes the form: 



TV 



N' 



-,(1x1) 



1=1 1=1 

N N' 

i<j li<h 
N 

i<j<s 



cr(^ + Ki + k'i) 



^(2X2) 



^^(0 ri/^l + lilf)(^{Kj + Kif) 



N' 



g»?!l+»)l2+')!3+Bii,!2+-E!i,l3+-Bi2.i3 
0-(C + + + Ks + + H'l, + '«i3) 



(3x3) 
hhhijs 



(2.69) 



where 



hh.-.lKijs. 



is the determinant of the N x N matrix of the selected entries of C. 

The dependence of r in (I2.69P on the discrete dynamical variables n, m, h comes only 
through the 9i and rij, which in turn depend on n, to, h via the plane wave factors (|2.7p and 



3 Reduction to the lattice KdV equations 

Dimensional reductions of the lattice KP systems can be obtained by imposing certain sym- 
metry conditions. This process could be done by imposing the condition that T^s ° Tgu — u 
for all lattice directions. This implies that the interchange 5 —S leads to a reversal of the 
lattice shift. 
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We now consider the plane-wave factors shifted in the " " direction 
We find that 



P 

V 



T~soTs{piy)^———^p:^. 



(3.1) 
(3.2) 

(3.3) 



If we take k = k' we have T^s ° Ts{pv) = pv. 
Let's consider the lattice KP 



w ~ w w — w w — w 



w — w w — w w — w 



(3.4) 



Setting A = —6, we have T\ o Tgw — w — w^^w — w^ where w — w{n — 1, to), yielding 



(w — w){w — w) ~ {w ~ w){w — w) 



(3.5) 



Similarly, setting A — —e, we have T\ o T^w = 'w = w^w~w^ where w = w{n,m — 1), 
yielding 



(w — w){w — w) = {{u — w){w — w) 



(3.6) 



This implies that the product on the right-hand sides of these relations is conserved in all 
directions, and then constant: 



{w — w){w — w) = constant , 
which is actually the elliptic version of the ^'lattice KdV" equation. 

From (I2.42al) . if we set X — —6 wc find that 

1 



(3.7) 



(^yj -w) = ^ { -paVa - Qa'Va 

Combining this equation with (I2.42bp and using the lattice KdV, we obtain 

PaiVcV^ - Vj^) = (?„(V„V„ - V„V„) , 
which is the "lattice modified KdV^^ equation. 



(3.8) 



(3.9) 



From (|2.53p . if we set A — —5, we get 

qsTT + psTT = psTT . 

and if we set A = — e, we get 

PSTT + qsTT = qsTT . 

p.lOp and (|3.1ip are the "lattice Hirota" equations. 



(3.10) 
(3.11) 



From (|2.62p . if we set A = —S we get 
and if we set A = —e we get 



(3.12) 



(3.13) 
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Using the identity: 



= , (3.14) 



we obtain the equation for Zp^a 

q^Zp^„ ~ q-pZp^^ p_pZp,^ - PaZfi^a 



(3.15) 



which is the ^Hattice Schwarzian KdV" equation. 

The dimensional reductions described here to lattice equations of the KdV class are 
obtained of the basis of the explicit form for the solutions described in the earlier sections. 
We would like to mention at this point that in the recent paper [2], it was shown how 
to embed all equations in the ABS list in the Schwarzian KP equation (jl.6p . through the 
canonical definition of a "Schwarzian variable" associated with each member of the list. 



4 The continuum limit 

In previous Sections, we derive the fully discrete KP equations. The derivation of interme- 
diate and full continuum limits of the lattice KP systems were studied in [221 [22] , as well 
as [3 [Sni 121] where the derivation of hierarchy arising from the lattice KP equations was 
systematically studied. We will consider continuum limit on the level of the elliptic soliton 
solution. We will give the detail only for the lattice KP (|2.49p and present the final results 
for the other lattice KP equations. 

We find out that it is more convenient to express the lattice KP equation (|2.49p in the 
following form 

(^0 - ^0 + C(e + e + A) - C(e + + A) - C(e) + CiS)) 
{wo -W0 + m + e) - C(e + C(e) + CiS)) 
(So - Wo + C(? + e + A) - C(e + ^ + A) - C(A) + C(5)) 

{wo -W0 + c(e + A) - c(e + s)- c(A) + as)) ' 

where wo = s^C(l + MC)^^r. 



(4.1) 



First continuum limit: The skew limit. We start to perform the contimmm limit 
by making change of discrete variables (n, m) — > {N = n + m, m) and then taking the limit 

e — (5 = ry— >0, m ^ oo ^ n— > — oo , s.t mr/ — t and N fixed . 

The plane wave functions (|2.7I) and (|2.8p become 

p{^) = (e-«^)>„«) " eC(«)-+C(^--)>,o,o(^^), (4.2) 

^{^') = (e^^^'^'pK') (e'^f'^"' e^(«)-'-«^+«>Po,o,o('*'), (4-3) 
and wqIti^ m, h] wo[N, r, h]. Then wq[N, t, h] satisfies the differential-difference equation 
# + p(e + ^ + A) - p{5) _ So - Wo + C(g + + A) ~ C(g + 2<5) - C(A) + m 



^ + p{(, + 5)- p{6) WO-WO + ai + A) - C(5 + 5)- C(A) + C(<5) 



(4.4) 



The second limit. Next, we perform the limit on the discrete variable h: 

A -> , ft, — > oo , s.t. J — Xh fixed. 
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The plane wave functions (j4.2l) and (14.31) become 

N 



Pin) - e«(«)''+«^-«)-';(«)>o,o.o(«), (4.5) 

K^') - {e^^'^'-' P.') e';«)'''-^('+''>+«''')>o,o,o('*'), (4.6) 

and wq [N, t, h] — > wq [N, t, 7] . Then wq [N^ t, 7] now satisfies the partial differential-difference 
equation 

- (1^ + p'io) = + p(0 " p(^)) (^^0 - 2z^o + !!^o + 2C(e) - c(e + 5) - c(e - <?)) • 

^ ^ ^ (4.7) 



The third Hmit: The Pull limit. This limit can be obtained by investigating the limiting 
behavior of the plane wave function. From (j4.5p , we have 

|^e-CW«p_^^ ^ gC(«)«+C('5-'^)T-CW7 

^ exp (iVln - (5) | - In |cr(K)| + C((5 - K)r - C(k)7 + C{£.)k) 
^ exp (^iV (^-5C('^) - yp('^) + y p'(«) + - C('«)7 

+r (^-c(aj) - 5p{K) ~ y p'(k) + .... j + c(e)«) 

^ exp (C(e)K + C(k)x + p(«;)y + p{n)'{K)t + ) , (4.8) 

where we define b — NS and 

X = — D — r — 7 , y = Td , and t = . (4.9 

/ . y 2 ' 6 2 ^ ' 

Applying the Taylor expansions, we have 

_ dwo d'^ d'^wo d^wo S'^ d'^wa 

wo = w, + S— + + y ^ + + ' ^4-10) 

where 6 = &o + -^(^ a-nd using the chain rule formulae, we have 

(4.ffa) 

Using (|4.10p and (|4.1ip . we recover the ''potential KP" equation in order 0((5^): 

3. , f 
2^0]^ + 2 

where wo[N,t,j] wo[x,y,t]. 



dwo _ 


dwQ 




dx 


dwo 


dwo 


dr 


dx 






djdr 


dx^ 



Wo]xt = 6[wo]2;[wo]2;i: + T;b"o]yy + T^i'^ohxxx , (4-12) 



At the end of this Section, we would like to mention that the continuum limit can be 
performed on the rest of KP equations namely the modified and Schwarzian KP equations. 
We will give a list of the results in the following 

The discrete modified KP equation: 



Vc 



, (4.13) 
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where Va = l—s'^C{l+MC) ^ (X-a,ic) ^ Taking the skew Umit, we have the differential- 
difference equation 



dr 



'-Vo. (xa,6Vo, - Xa,xVc}j + (xa.A^^ - Xcc,sV o) = , (4.14) 

where Va[n,m,h] — >■ Va [N, r, h] . 



Taking another hmit on the discrete variable h, we obtain the partial differential-difference 
equation 

dldT Va ^ Va^ ^ 

where Xa = C(a) + C(C + ^) - C(C + " + 'J) and y„ [TV, r, /i] ^ F„ [AT, r, 7] . 
Taking the full limit, we obtain the modified KP cqiiation 

[VaU = [VaU. + W<.]vv - ^KUVaU - Q[Va\y[VaU , (4.16) 

where V^ [N, r, 7] V^ [x, y,t]. 

The lattice Schwarzian KP equation: 

(1 - Xa,sSl3,a - Xp-sSfi,o^ (l " Xa,A^/3,a - %,-\Sl3,c^ - Xa,eSp,a - X/S-e-^/J.a^ 

(1 - Xa,\S/3^a - XlS-xSi3,a^ ^1 - Xa,6Si3,a " X/3,-5'5/3,a^ (l " Xa,eSl5,a " X/3,-e'S'/3,a) 

= 1 , (4.17) 

where Sp,a = S"^ [x,3,a:']~^ C [1 -|- MC]~^ [Xa,K]~^ r. Taking the skew limit, we derive the 
differential-difference equation 

ipil + a) - p{0)S^,a + ipj^+P) - P(|))^^,a +Xa,A^ +f/3,-A%^ 
1 - Xa,\Si3^a - Xfi,-\Sp,a 

1 — Xa,sS/3,a - Xl3,-S^0,a 

where S^^ain, m, h] — > Sjs^aW, t, h]. 

Taking another limit on the discrete variable h, we have the partial differential-difference 
equation 

[P(C + /3) + P(e + a) - 2p(0]5;3,„ + [Xa,* + X^,-5] 



9r 9r97 



1 - [Xa,5 + Xp-5]Sff,a + 



07 

dSa 



^ [p{6) - p( 0]5ff,a + [p((5) + p($ + /3)]5;3,„ + Xa,5^ 



1 - Xa,5Sp^a — XP^-sSp^a 

^ [p{i + a)- p{.5)]Sf3,a - [p{5) + p{i)]SfS,a + Xff,-^-9:r ^ q (-4 j_g) 

1 - Xa.^S/S.a - X/3,-i5<S'/3,a 



9S/3,o 
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where Sp^a [N, r, /i] Sp^a [N, r, 7] . 



Taking the fuU hmit, we obtain the Schwarzian KP equation 

[Sli^a\t ~ [Sp,a\xxx '^[^ [i ,a\'y ~ ^[^ P ,a\^x\ 9 ( [S p ^a\y 



where Sp^a[N,T,h] Sp^a[x,y,t]. 



= , (4.20) 



5 Summary 

In this paper, we estabhshed the expUcit form of a class of eUiptic sohton solutions for all 
the lattice KP equations, based on a construction using elliptic Cauchy matrices denoted 
by M. urthermore, the construction exhibits numerous relations between the various lat- 
tice equations, as well as corresponding Lax pairs. The explicit form for the corresponding 
r-function depends crucially on the coefficient matrix C, which opens the way to classify 
the various different lattice soliton behaviours according to the Schubert decompositions 
of the corresponding Grassmannians, following similar work in the continuous case by Ko- 
dama and Chair avarty, |16| . Several reductions were considered: i) dimensional reduction 
to KdV lattice systems, and ii) continuum limits to the semidiscrete and fully continuous 
KP equations. For all these equations the corresponding elliptic soliton solutions are derived 
in parallel. Furthermore, the result in (|2.69p can be simplified to the cases of trigonomet- 
ric/hyperbolic by taking: a{x) t-^ sin(a;) or <7{x) n> sinh(x) and for the rational case we have 
a{x) H» X. 

There are well-established connections between soliton solution of integrable PDF and 
the integrable many body systems [261 127| . This can be made most explicit in he rational 
and trigonometric/hyperbolic cases. The elliptic case of this correspondence is more difficult 
to establish, but we expect that the elliptic solitons which we have studied here can be 
connected to elliptic case of the discrete-time Ruijsenaars model constructed in [53], but at 
this juncture this is still conjectural. Nevertheless, in [28l[29|, an explicit connection between 
the rational discrete-time RS system and the KP lattice was established using the solution 
structure of the former model. 
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